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摘 m. 研究 大 气 水 汽 的 反 演 及 了 解 其 时 空 变化 ,对 有 效 讨 
地 区 CE -318 观测 的 大 气 水 汽 数据 ,结合 探 空 水 汽 数据 和 地 耳 


F 估 大 气 生态 服务 功能 具有 重 


现实 意义 。 基 于 石家庄 
水 汽 压 数 据 ,运用 传统 回归 分 析 、 改 进 型 回归 分 析 、 


分 段 反 演 等 方法 ,构建 传统 回归 模型 .基于 CE - 318 的 大 气 水 汽 改进 型 模型 和 分 季节 模型 3 种 大 气 水 汽 反 演 模型 ， 


并 经 过 精度 对 比 评估 选 出 适用 于 研究 


的 大 气 水 汽 最 优 反 演 模型 。 检 验 结果 表明 ,基于 CE -318 的 大 气 水 汽 改进 
个 


模型 中 的 最 佳 值 ,是 该 地 区 最 优 的 大 气 水 汽 反 演 模型 。 


[X 
型 反 演 模型 的 各 项 精度 检验 指标 均 为 3 
关键 词 : 大 气 水 汽 ; CE -318; 探 空 ; 地 


大 气 水 汽 含量 是 进行 遥感 定量 反 演 工作 必须 要 
考虑 的 一 个 关键 参数 中 ,同时 水 汽 循 环 又 是 预测 气 
候 变 化 的 一 个 重要 参量 ,影响 着 地 表 气 候 的 演变 , 因 
此 ,获取 连续 性 的 大 气 水 汽 数据 始终 是 相关 领域 探 
讨 的 热门 话题 。 目 前 ,在 技术 方式 上 ,大 气 水 汽 的 反 
演 方法 分 为 遥感 法 和 实测 法 ,遥感 法 又 可 根据 设备 
安放 位 置 的 不 同 可 分 为 卫星 遥感 和 地 基 遥 感 。 卫 星 
遥感 的 优势 是 可 获取 大 面积 的 水 汽 时 空 分 布 信息 ， 
且 卫 星 遥 感 水 汽 产品 众多 号 ;目前 学 者 更 倾向 于 多 
源 遥 感 数据 一 体 化 来 反 演 水 汽 , 如 王 永 前 等 2 基于 
物理 统计 思想 设计 了 一 个 融合 多 源 遥 感 数据 的 水 汽 
反 演 算法 ,该 算法 极 大 地 缩减 了 多 云天 气 条 件 的 影 
响 ; 此 外 ,对 卫星 降水 产品 的 校正 中 也 是 该 方向 的 
一 个 热门 研究 话题 。 地 基 遥 感 主要 有 地 基 GPSS A 
地 基 CE -318( 太 阳光 度 计 ) 等 ,运用 这 些 方法 的 优 
势 在 于 受 地 表 辐 射 和 反射 的 影响 很 小 ,探测 精度 更 
高 ;其 中 地 基 CE -318 太阳 仪 是 近 几 年 地 基 遥 感 中 
的 一 种 简便 的 常规 探测 方法 。 利 用 其 936 nm 水 
汽 吸收 通道 所 获 太 阳 辐 射 数据 ,基于 改进 型 Langley 
算法 反 演 的 水 汽 具 有 监测 精度 高 .时 间 分 辩 率 高 
和 耗 能 低 的 优势 。 实 测 法 主要 是 高 空 资料 计算 ， 
即 利用 探 空气 球 测 出 大 气 中 每 一 层 的 湿度 .气温 、 气 
压 等 气象 因素 ,采用 空间 外 推 和 内 插 算 法 推算 出 大 
气 水 汽 的 含量 。 运 用 这 些 方式 研究 大 气 水 汽 的 反 
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演 , 了 解 其 时 空 变化 规律 ,对 监测 周边 大 气 水 汽 环境 
演变 .加 强生 态 环境 建设 具有 重要 的 现实 意义 。 

d MESE HH UT ,大 气 水 汽 与 地 面 水 汽 压 二 者 
拥有 很 好 的 关联 性 ,可 利用 两 者 之 间 相 应 的 规律 性 
建立 关系 模型 ,成 为 一 种 新 的 估算 大 气 水 汽 的 方法 ， 
估算 精度 也 很 高 。 王 研 峰 等 "探讨 了 高 原 地 区 上 
空 水 汽 含量 与 地 面 水 汽 压 之 间 的 相关 性 ,发 现 二 者 
之 间 存在 二 次 多 项 关系 ,并 证 实 该 关系 式 具 备 一 定 
的 实用 价值 。 姚 俊 强 等 中 -中 以 天 山 山区 为 研究 区 ， 
运用 探 空 水 汽 数据 建立 了 与 地 面 水 汽 压 的 关系 式 ， 
并 推算 了 天 山 山区 的 水 汽 含量 ,解析 了 该 地 域 大 气 
水 汽 的 时 空 布局 及 其 与 气象 要 素 的 关系 。 潘 硝 
等 ' 引 基于 统计 思维 探究 了 大 气 水 汽 与 地 面 水 汽 压 之 
间 的 相关 性 ,构建 出 二 者 关系 计算 模型 ,并 以 GIS 为 
分 析 平 台 , 完 成 了 对 大 气 水 汽 空 间 布 局 的 仿真 模拟 。 
张 丹 等 “中 利用 MODIS 反 演 的 水 汽 数 据 和 地 面 常规 
气象 材料 ,构建 了 中 国 以 逐 月 地 面 水 汽 压 来 估算 大 气 
水 汽 的 统计 模型 ,并 对 模型 精准 度 做 了 相应 误差 检 
验 ,结果 证 明 该 模型 算法 是 一 种 切实 可 行 的 方法 。 其 
他 学 者 5 -9 也 在 相关 领域 做 了 深入 的 研究 。 

鉴于 前 人 所 做 地 面 水 汽 压 反 演 大 气 水 汽 模 型 的 
研究 工作 ,可 发 现 其 大 部 分 采用 传统 经 验 模型 法 进 
行 反 演 分 析 , 这 使 模型 反 演 结果 存在 一 定 程度 的 精 
度 误 差 。 另 外 ,模型 反 演 精度 也 会 受 地 理 环境 复杂 
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I) 


于 时 


程度 和 和 气候 状况 的 影响 较 大 ,因此 仍 需 根据 研究 区 
的 实际 情况 ,构建 适合 研究 区 的 反 演 模型 。 本 文通 
过 改进 型 回归 模型 法 ,结合 2016 年 10 月 至 2017 年 
9 月 石家庄 地 区 的 CE -318 观测 大 气 水 汽 数据 、 探 
空 水 汽 数 据 和 地 面 水 汽 压 数据 ,构建 石家庄 地 区 上 
空 大 气 水 汽 反 演 模型 ,以 提供 适用 于 该 地 区 估 测 大 
气 水 汽 含量 的 新 方法 ,为 监测 空中 大 气 水 汽 含量 的 
长 时 间 序 列 变化 提供 依据 。 


1 数据 来 源 及 处 理 


数据 来 源 

采用 的 CE -318 观测 数据 源 于 可 自动 工作 的 
新 型 CE - 318 地 基 太 阳光 度 计 所 获取 的 数据 ,该 
仪 需 在 936 nm 波段 处 有 吸收 水 汽 的 特性 ,联合 
870 nm 和 1 020 nm 波段 通道 可 反 演 出 大 气 柱 水 汽 
含量 5 。 仪 器 架设 点 地 理 坐标 为 38.00*N， 
114.52"E ,经 定 标 和 安装 后 正常 工作 ,观测 时 间 为 
日 升 到 日 没 , 北 京 时 冬 半 年 为 08 :00 一 18 :00 , 夏 半 
年 7:00 一 19:00。 观 测 所 获 数据 均 为 1.5 级 去 云 高 
质量 数据 。 

探 空 数据 来 自 美国 NCDC (The National Climat- 
ic Data Center) 的 IGRA Zt fis JÆ ( ftp ://ftp. nede. no- 
aa. gov/pub/ data/ghcn/daily/ (2017 - 10 - 20) J ,该 
数据 集 已 经 过 严格 质量 控制 。 考 虑 到 石家庄 地 
区 无 探 空 站 点 ,经 对 比 地 形 地 貌 和 气候 等 自然 条 件 
相对 一 致 的 情况 下 ,从 数据 集中 下 载 距 离 石 家 庄 水 
平 距离 最 近 和 海拔 高 度 相 差 无 几 的 位 于 37. 18? N, 
114.35? E. 的 邢台 探 空 站 所 获 资 料 进行 水 汽 反 演 。 
探 空 仪器 观测 时 间 只 有 UTC 0 时 和 UTC 12 时 , 数 
据 包含 的 变量 有 各 大 气压 层 的 气温 \ 温 度 露 点 差 HC 
压 、 风 速 和 风向 等 ,根据 这 些 变量 利用 探 空 积分 法 可 
反 演 计算 出 整个 垂直 大 气 层 的 大 气 水 汽 值 "” 。 

地 面 水 汽 压 数据 源 于 国外 Meteomanz. com 服务 
网 ,该 网 站 可 提供 来 自 全 球 各 气象 站 点 官方 发 布 的 
气象 数据 ,从 中 下 载 石家庄 地 区 2016 年 10 月 到 
2017 年 9 月 的 小 时 级 (UTC 00, UTC 03 、UTC 06, 
UTC 09 , UTC 12 .UTC 18 和 UTC 21) 变量 数据 ;之 后 
利用 马 格 纳 斯 经 验 公 式 , 根 据 所 获 气温 和 相对 湿度 
两 种 气象 要 素 计 算出 当地 水 汽 压 值 。 考 虑 到 数据 的 
准确 性 ,笔者 将 中 国 气象 数据 网 中 公布 的 2000 年 1 
月 至 2017 年 3 月 月 均 水 汽 压 值 与 Meteomanz. com 
网 站 上 相应 计算 出 的 月 均 水 汽 压 值 ( 共 207 组 数 


1.1 


HE) 进行 了 对 比分 析 , 发 现 它们 的 相关 系数 为 
0.999 7, 在 0.01 的 置信 区 间 上 极 显著 相关 ,并且 通 
过 了 下 和 了 检验 ,说 明 该 网 站 上 数据 准确 度 和 可 靠 
性 很 高 ,可 用 来 做 科学 研究 。 
1.2 数据 筛选 处 理 

经 过 晴好 少 云 天 气 下 数据 量 的 筛选 ,研究 期 间 
2016 年 10 月 至 2017 年 9 月 CE -318 观测 的 大 气 
水 汽 可 用 数据 量 为 147 d。 但 是 由 于 CE -318 观测 
受 环境 限制 ,在 落日 与 阴雨 天 气 或 空气 质量 很 差 时 ， 
仪器 停止 观测 ;并 且 UTC 12 时 所 获 探 空 数 据 等 于 
石家庄 地 区 北京 时 20 时 的 数据 ,太阳 此 刻 即 将 或 早 
已 落 山 ,CE -318 仪器 暂停 工作 。 此 外 ,本 文 获取 的 
地 面 水 汽 压 同样 是 以 世界 时 为 基准 进行 计算 的 。 考 
虑 到 3 种 数据 在 研究 时 段 一 致 性 的 问题 ,只 选取 
UTC 0 时 刻 数 据 作 为 研究 样本 。 在 进行 UTC 0 时 
CE -318 观测 水 汽 转 化 和 计算 时 ,依照 探 空 探 测 时 
间 大 约 1p 为 准 ,选取 CE -318 仪器 在 当天 23 :30 
至 次 日 0:30 反 演 的 大 气 水 汽 取 均 值 , 得 到 UTC 0 时 
的 CE -318 实测 水 汽 值 。 


2 大 气 水 汽 反 演 模型 的 构建 


经 过 数据 异常 值 的 吻 除 和 筛选 , 共 提取 出 127 
组 同一 时 刻 的 数据 用 于 模型 构建 。 同 时 ,将 127 组 
数据 分 为 2 组 :构建 模型 数据 和 模型 检验 数据 ;分 组 
依据 为 :根据 CE -318 观测 水 汽 值 将 全 部 数值 重新 
排序 ,在 排列 数值 中 从 低 到 高 每 隔 4 个 数值 抽取 1 
个 数值 用 来 误差 检验 ,这 样 原 始 数据 被 分 成 102 组 
建 模 数 据 及 25 组 模型 检验 数据 。 

2.1 传统 回归 模型 法 构建 CE -318 大 气 水 汽 反 演 
模型 

以 建 模 数据 中 地 面 水 汽 压 为 模型 中 的 自 变量 ， 
CE - 318 观测 水 汽 为 模型 中 的 因 变 量 ,用 统计 软件 
建立 二 者 的 6 种 估算 模型 类 型 ,分 别 为 指数 函数 线 
性 函数 .对 数 函 数 . 乘 过 函数 二 次 函数 和 三 次 函数 
( 表 1)。 

从 表 1 可 以 看 出 ,6 种 回归 模拟 方程 都 通过 了 
0.01 的 显著 性 FF 检验。 相 比 决定 系数 (R ) 可 知 ,二 
次 函数 和 三 次 函数 的 R 最 高 且 相同 , 均 为 0. 924 2, 
本 着 计算 简便 性 的 原则 ,选取 二 次 函数 作为 地 面 水 
汽 压 (Eng) 与 CE -318 观测 水 汽 (Pu -s1 ) 的 最 
优 关系 模型 选择 ,表达 式 为 : 

PW, 44 =0.001 6E, +0. 101 4E, +0.280 3 (1) 


武 紫 超 等 :石家庄 地 
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区 大 气 水 汽 的 反 演 模型 


表 1 CE -318 观测 水 汽 与 地 面 水 汽 压 回归 方程 拟 合 结果 
Tab.1 The fitting results of water vapor and surface vapor pressure observed by CE -318 


模型 类 型 可 归 模 拟 方程 R? r SE F a =0.01 
指数 函数 y 20. 476 4e0.089 2 0.855 8 0.925 1 0. 048 3 541.99 P «0.01 
线性 函数 y 20. 146 4x +0.103 5 0.917 4 0.957 8 0.0354 1 110. 20 P «0.01 
对 数 函 数 y=1.221ln(x) -0. 801 4 0. 769 6 0. 877 3 0. 059 0 334. 05 P «0.01 
TETE PRU y 20. 2349 8551 0. 904 2 0. 950 9 0.037 5 975.75 P «0.01 
次 函数 y 20. 001 6a? +0. 101 4x +0. 280 3 0. 924 2 0. 961 4 0. 033 9 1 220. 06 P «0.01 
三 次 函数 y= -1x107532 +0.0017xz2+0.1008x+0.2818 0.924 2 0. 961 4 0. 033 8 1 220. 12 P «0.01 
注 :y 为 Per 318 值 ,x 为 Eug E; SE 为 标准 误 。 下 同 。 
AA ， 探 空 观测 水 汽 与 地 面 水 汽 压 回归 方程 拟 合 结果 
Tab.2 The fitting results of water vapor and surface vapor pressure from sounding observation 
模型 类 型 可 归 模 拟 方程 R? r SE F a=0.01 
IEE y 20.354 8e” 1% 6: 0. 846 3 0. 919 9 0. 061 0 421.38 P «0.01 
线性 函数 y 20. 165 8x -0.101 0.917 5 0.957 9 0. 040 0 1 112.13 P «0.01 
对 数 函 数 y=1.378 3In(x) -1.1172 0.764 8 0. 874 5 0. 067 5 325.19 P «0.01 
乘 军 函数 y=0. 143 4« $3 0.927 6 0.963 1 0. 039 8 1 124.25 P «0.01 
次 函数 y 20. 002x? +0. 110 5x +0. 116 3 0.925 6 0. 962 1 0. 038 0 1 244.35 P «0.01 
三 次 函数 y= -3x10-5x3 +0.002 2x? +0. 108 9x +0. 120 5 0. 925 6 0. 962 1 0. 038 0 1 243. 76 P «0.01 
注 :y H PW. ff, x 为 Bw 面值 。 
2.2 ”改进 型 回归 模型 法 构建 CE -318 KAKAR ”其 表达 式 如 下 : 
1.0163 
演 模型 PW... =0. 143 AES (2) 


2.2.1 基于 探 空 观测 水 汽 的 常规 反 演 模型 ”建立 
探 空 观测 水 汽 与 地 面 水 汽 压 之 间 相 关联 程度 ,构建 
基于 探 空 观测 水 汽 的 常规 反 演 模型 。 运 用 统计 回归 
分 析 方 法 ,在 置信 和 度 为 99% 的 区 间 内 ,利用 线性 函 
数 .指数 函数 . 乘 寡 函数 .对 数 函 数 .二 次 函数 和 三 次 
函数 等 6 种 方程 类 型 , 拟 合 探 空 观测 水 汽 (PWjs ) 
与 地 面 水 汽 压 (Ej ) 的 回归 表达 式 ( 表 2)。 


2.2.2 基于 探 空 与 CE -318 的 大 气 水 汽 修正 模型 

在 进行 修正 模型 构建 之 前 需 对 CE - 318 和 探 空 
两 种 探测 手段 的 水 汽 计算 结果 进行 对 比分 析 , 以 确 
保修 正 模型 建立 的 科学 性 和 严谨 性 。 图 1 分 别 是 两 
种 观测 手段 计算 水 汽 结果 的 趋势 图 和 拟 合 图 。 由 图 
la 明显 看 出 ,二 者 大 气 水 汽 值 随 时 间 演 变 趋势 基本 
一 致 ,峰值 和 谷 值 接近 ,说明 两 种 手段 均 能 观测 出 大 


Hx 2 结果 可 看 出 , 探 空 观测 水 汽 与 地 面 水 汽 
压 二 者 具有 很 好 的 相关 性 , 且 6 个 方程 的 拟 合 结 
都 能 通过 0. 01 极 显 著 水 平王 检验 ,符合 统计 学 要 
ok. Hen R 为 0.927 6 的 乘 窜 函数 拟 合 效果 最 好 ， 
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气 水 汽 的 时 间 变 化 信号 。 不 过 ,二 者 大 气 水 汽 在 具体 
数值 上 有 一 定 的 差异 ,系统 误差 和 偶然 误差 的 存在 以 
及 两 种 观测 仪器 所 处 位 置 不 同 , 相 对 地 势 和 海拔 
有 一 定 差距 ;此 外 ,同一 时 刻 两 种 仪器 所 在 地 气象 条 
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注 :PWCE a8 为 CE - 318 观测 水 汽 , PW pes 为 探 空 观测 水 汽 。 下 同 。 
图 1 CE -318 与 探 空 计算 所 得 水 汽 值 趋势 图 (a) 和 散 点 图 (b) 


Fig. 1 


Compared results of water vapor charts and scatter charts between CE —318 and sounding observations 


干 a 


RI 探 空 观测 水 汽 与 CE -318 观测 水 汽 回 归 方 程 拟 合 结果 
Tab.3 The fitting results of water vapor from sounding and CE —318 observations 


模型 类 型 可 归 模 拟 方程 R r SE F a 20.01 
EEE y 20. 526 740-508 6: 0. 833 2 0. 912 8 0. 061 5 300. 56 P «0.01 
线性 函数 y 20. 862 3x +0. 225 1 0. 953 4 0. 976 4 0.0265 2047.07 P «0.01 
对 数 函 数 y=1.194 9In(x) 41.531 7 0. 820 8 0. 906 0 0.052 1 458. 08 P «0.01 
RE BAUR y 21. 13549 8065 0. 939 0 0. 969 0 0. 027 0 1 980. 27 P «0.01 
次 函数 y= -0.015 3x? +0.936 4x +0. 176 9 0.954 1 0. 976 8 0.0263 2079.74 P «0.01 
次 函数 y 20. 000 2x3 -0.016 6x? +0.939 2x +0. 175 7 0.954 1 0. 976 8 0.0063 — 2079.72 P «0.01 
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件 不 同等 ,都 会 导致 误差 的 出 现 。 由 图 1b 显示 ,二 
者 紧密 分 散在 一 条 直线 周围 ,R 达到 了 0.952 3 ,可 
见 两 者 相关 度 很 高 ,说 明 它 们 之 间 拥 有 很 好 的 线性 
关联 性 。 

同样 ,基于 统计 回归 分 析 法 对 探 空 与 CE -318 
观测 水 汽 2 组 数据 进行 6 种 回归 模型 的 拟 合 , 拟 合 
情况 见 表 3。 

从 表 3 中 RT 能 够 看 出 ,前 4 种 模型 对 两 者 的 拟 
合 效 果 都 不 太 理 想 , 而 R 为 0.954 1 的 二 次 函数 和 
三 次 函数 拟 合 度 最 高 ,而 且 都 通过 了 0.01 极 显著 水 
平 的 了 检验 ,同样 出 于 计算 简便 性 考虑 ,选取 二 次 
函数 方程 作为 PW 探 空 与 CE - 318 观测 水 汽 的 回 
归 模 型 ,其 表达 式 为 : 

PW as = -0.015 3PW; 


4 +0. 936 4PW gx + 


空 RE 

0. 176 9 (3) 
2.2.3 CE -318 大 气 水 汽 反 演 模型 的 建立 ”将 式 
(2) 和 式 (3 ) 进 行 整合 ,得 到 地 面 水 汽 压 与 CE -318 
观测 水 汽 的 相关 关系 ,构建 CE - 318 大 气 水 汽 反 演 
模型 ,从 而 可 以 对 长 时 间 序 列 的 大 气 水 汽 动态 进行 
监测 ,为 合理 利用 大 气 水 汽 提供 决策 支持 。 地 面 水 
汽 压 与 CE -318 观测 水 汽 的 回归 模型 表达 式 为 : 


PW 4.2 -0.000 3E © +0. 134 3Egu + 


0. 176 9 (4) 
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对 地 面 水 汽 压 按 模 型 进行 计算 ,得 到 CE - 318 
大 气 水 汽 模 拟 值 ,将 CE - 318 大 气 水 汽 模 拟 值 与 
CE -318 大 气 水 汽 实 测 值 进行 对 比 ,结果 见 表 5。 各 


通过 该 回归 模型 计算 得 出 的 CE -318 大 气 水 
汽 模 拟 值 ,与 CE -318 大 气 水 汽 实测 值 进行 对 比 ， 
二 者 R 达到 了 0.914 8 ,标准 误差 0.364 3 , 拟 合 效 
果 较 为 理想 。 
2.3 ”分 季节 构建 CE -318 大 气 水 汽 反 演 模型 

考虑 到 由 于 运用 单一 长 时 间 序 列 地 面 水 汽 压 估 
算 大 气 水 汽 含量 的 精准 度 问 题 , 通 过 传统 四 季 划 分 
的 指标 ,将 CE -318 大 气 水 汽 取 值 划分 为 春 、 夏 、 
秋冬 4 个 区 间 段 ( 表 4) ,在 不 同 的 区 间 段 内 选取 相 
应 的 地 面 水 汽 压 ,按照 上 述 改进 型 回归 模型 法 的 原 
理 进行 CE -318 大 气 水 汽 反 演 模型 的 构建 。 


表 4 分 季节 模型 下 建 模 数据 和 检验 数据 统计 量 
Tab.4 Statistic results of modeling data and test data 


under the seasonal models 


分 段 类 别 建 模 数据 量 模型 检验 数据 量 
春 34 9 
夏 22 5 
$k 23 5 
A 23 6 
总 计 102 25 


在 构建 模型 时 ,同样 在 不 同 季节 区 间 段 内 分 别 
按照 指数 函数 .线性 函数 .对 数 函 数 . 乘 寡 函数 二 次 
函数 和 三 次 函数 6 种 方式 进行 模型 的 构建 ,之 后 根 
JE R 选取 各 季节 区 间 段 内 的 最 优 模型 为 : 


( 春 ) 
(E) (5) 
( 秋 ) 
( 冬 ) 


[ELS R^ 达到 了 0. 93 ,标准 误差 为 0.330 2。 
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表 5 分 季节 区 间 回 归 模 型 法 拟 合 结果 


Tab.5 The fitting results of seasonal interval regression model 


季节 区 间 最 优 模型 R? r SE F a 20.01 
春 y=1.012 1x 40. 002 3 0. 856 1 0.925 3 0.038 5 190. 42 P «0.01 
夏 y=1.140 3x -0.5649 0.702 8 0. 838 4 0.1120 47. 30 P «0.01 
秋 y=1. 128 2x -0.1326 0. 804 9 0. 897 2 0. 069 2 86. 64 P «0.01 
A y 21.065 7x —0.017 5 0. 583 1 0. 763 6 0. 032 2 29. 38 P «0.01 
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3.1. 建 模 前 后 大 气 水 汽 与 地 面 水 汽 压 的 相关 性 比较 

通过 比较 建 模 前 和 建立 不 同 模 型 后 的 大 气 水 汽 
与 地 面 水 汽 压 之 间 的 相关 性 ,可 直观 地 观察 出 各 模 
型 的 优 劣 情况 ,因此 ,从 比较 建 模 前 大 气 水 汽 与 地 面 
水 汽 压 数 据 、 建 立 传统 模型 及 改进 模型 后 的 大 气 水 
汽 与 地 面 水 汽 压 数据 .划分 季节 后 的 大 气 水 汽 与 地 
面 水 汽 压 数 据 之 间 的 相关 性 来 研究 大 气 水 汽 与 地 面 
水 汽 压 之 间 的 关系 ,结果 如 表 6 所 示 。 


表 6 石家庄 地 区 建 模 前 后 大 气 水 汽 与 地 面 水 汽 压 
相关 性 对 比 

Tab.6 Compared results of the correlation between 

atmospheric water vapor and surface vapor pressure 


before and after modeling in Shijiazhuang area 


e 建 模 前 建 模 后 PW - Ey 
s PW -Eww ”传统 模型 。” 改进 模型 
全 年 0.957 8 0.996 3 0.999 8 
春季 0.922 3 0.994 7 0.998 6 
夏季 0.838 6 0.979 8 0.990 0 
秋季 0.905 7 0.980 3 0.998 0 
冬季 0.761 8 0.999 7 0.999 5 


注 :P 多 为 大 气 水 汽 ,5 而 为 地 面 水 汽 压 。 


由 表 6 可 得 ,全 年 及 不 同 季节 建 模 后 大 气 水 汽 
与 地 面 水 汽 压 的 相关 性 均 有 所 提高 ,其 改进 模型 所 


| y=1.083 6x-0.0269 oo o 


得 大 气 水 汽 与 地 面 水 汽 压 的 相关 性 除 冬季 外 均 为 最 
高 ,全 年 两 者 的 相关 值 无 限 接近 于 1, 可见 ,改进 模 
型 更 能 体现 大 气 水 汽 与 地 面 水 汽 压 之 间 的 关系 。 
3.2 ”模型 精度 检验 

为 了 探讨 各 模型 反 演 大 气 水 汽 的 情况 , 需 对 其 
进行 精确 度 评估 ,以 检验 模型 的 有 效 性 。 用 于 检验 
的 数据 包括 地 面 水 汽 压 值 和 CE - 318 实测 水 汽 值 。 
将 地 面 水 汽 压 数据 分 别 带 入 到 各 模型 中 进行 计算 ， 
得 到 CE -318 大 气 水 汽 模拟 值 ,进而 可 利用 CE - 
318 大 气 水 汽 模拟 值 与 实测 值 的 相关 关系 来 检验 和 
评价 模型 的 反 演 精度 。 根 据 6 个 指标 , 即 总 体 预测 
精度 R 、 标 准 误 (SE) \ 均 方 根 误差 (RMSE) 平均 标 
准 误差 (MSE) 和 平均 绝对 误差 (MAE) 对 模型 精度 
进行 对 比 说 明 。 各 模型 拟 合 结果 及 精度 对 比 结果 如 
图 2 和 表 7 所 示 。 

结果 显而易见 ,3 个 模型 中 实测 值 与 模拟 值 的 
相关 程度 较 好 , 尼 都 在 0.8 以 上 ,由 此 可 知 ,这 3 个 
模型 的 预测 精准 度 都 较 高 ,具有 一 定 的 可 行 性 和 科 
学 性 。 通 过 3 个 模型 的 各 个 指标 对 比 可 发 现 , 在 总 
体 预测 精度 比较 上 ,改进 型 回归 模型 的 精准 度 最 高 ， 
其 次 为 分 季节 回归 模型 和 传统 回归 模型 ;在 R^ SE 
fI RMSE 比较 上 ,改进 型 回归 模型 的 指标 值 最 低 ,其 
次 为 传统 回归 模型 和 分 季节 回归 模型 ;在 MSE 和 
MAE 的 指标 对 比 上 ,同样 是 改进 型 回归 模型 的 指标 
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图 2 不 同 模型 下 CE -318 大 气 水 汽 实测 值 与 模拟 值 相关 性 分 析 


Fig. 2 


Analyzed results of the correlation between measured and simulated values of atmospheric water vapor observed by CE -318 


under different models 


R7 各 模型 精度 对 比 


Tab.7 Comparison of the model accuracy 


模型 总 体 预测 精度 /% R? MSE SE RMSE MAE 
传统 回归 模型 76.72 0.883 0 0.141 3 0.082 0 0.409 8 0.282 5 
改进 型 回归 模型 79.06 0.892 7 0.137 5 0.078 5 0.405 0 0.275 0 
分 季节 回归 模型 77.65 0.822 5 0.153 8 0.101 0 0.5107 0.307 6 


值 最 低 ,其 次 为 传统 回归 模型 和 分 季节 回归 模型 。 
所 以 ,通过 各 指标 对 比 得 出 结果 为 :改进 型 回归 模型 
的 各 指标 均 最 佳 ,精度 最 高 ,可 以 作为 石家庄 地 区 上 
空 大 气 水 汽 反 演 的 最 优 估算 模型 。 
3.3 ”模型 适用 性 分 析 

以 北京 地 区 39. 992°N,116. 3109 E 为 例 ,对 该 
模型 进行 普 适 性 分 析 。 所 用 CE -318 水 汽 数据 源 
于 美国 AERONET 全 球 气 溶胶 探测 网 站 ,地 面 水 汽 
压 数 据 源 自 Meteomanz. com 服务 网 。 同 样 ,按照 时 
间 一 致 性 原则 经 过 数据 筛选 处 理 之 后 ,将 UTC 0 时 
的 地 面 水 汽 压 值 代入 改进 型 模型 中 ,得 出 的 水 汽 结 
果 与 实测 结果 进行 对 比分 析 , 用 来 检验 该 模型 的 普 
适 性 。 对 比 效果 如 表 8 和 图 3 所 示 : 


WŒ MSE 为 平均 标准 误差 ,SE 为 标准 误 ,RMSE 为 均 方 根 误差 ,MAE 为 平均 绝对 误差 。 
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图 3 CE -318 大 气 水 汽 实 测 值 与 模拟 值 
相关 性 分 析 


Fig.3 Analyzed results of correlation between the measured 
and simulated values of atmospheric water vapor from 


CE - 318 observation 


表 8 模型 相关 性 及 精度 对 比 


Tab.8 Comparison of model correlation and accuracy 


相关 性 (全 年 ) 对 比 结果 精度 对 比 
建 模 前 PW -Epam 建 模 后 PW-Ewg 总 体 预测 精度 /% R? SE RMSE 
0.955 4 0.991 7 传统 模型 70.74 0.952 5 0.032 5 0.384 8 
0.999 7 改进 模型 74.11 0.959 0 0.030 2 0.589 2 
在 相关 性 比较 上 ,北京 地 区 建 模 前 后 大 气 水 汽 。 子 , 长 时 间 监 测 其 动态 变化 ,对 生态 系统 服务 评价 、 


与 地 面 水 汽 压 相关 性 呈现 效果 与 石家庄 地 区 一 致 ， 
都 是 改进 型 模型 两 者 的 相关 性 最 高 , 极 接近 于 1。 
其 次 ,从 两 种 水 汽 值 的 精度 对 比 和 拟 合 效果 上 也 可 
看 出 , 除 均 方 根 误差 指标 上 改进 模型 较 劣 于 传统 模 
型 之 外 ,其 他 指标 均 优 于 传统 模型 ,尤其 是 总 体 预 测 
精度 高 于 传统 模型 近 4 个 百分点 ,这 说 明 该 模型 也 
可 很 好 地 适用 于 北京 地 区 。 人 研究 效果 表明 该 模型 可 
推广 到 其 他 地 区 ( 如 京 津 旨 地 区 ) 使 用 ,具有 一 定 的 
普 适 性 。 但 受 地 理 环境 和 气候 因素 的 影响 ,在 使 用 
时 需 根据 地 区 实际 情况 ,对 模型 参数 稍 加 修改 , 便 可 


适用 于 所 有 有 需求 的 地 区 。 
4 ”结论 与 讨论 


环境 保护 和 社会 生产 发 展 具 有 重要 的 参考 价值 。 
笔者 首先 建立 探 空 水 汽 与 地 面 水 汽 压 的 常规 水 汽 
模型 ,再 利用 CE -318 水 汽 与 探 空 水 汽 之 间 的 关系 
进行 修订 ,从 而 构建 基于 CE -318 的 大 气 水 汽 改进 
型 反 演 模型 ,并 进一步 依照 四 季 构 建 分 季节 的 大 
气 水 汽 反 演 模型 ,最 终 运用 相关 指标 对 模型 做 精准 
度 检 验 对 比 , 选 出 适合 该 地 区 的 最 优 反 演 模型 。 
利用 研究 数据 构建 出 传统 回归 模型 [公式 
(1)]、 基 于 CE -318 的 大 气 水 汽 改进 型 反 演 模型 
【公式 (4)] 及 分 季节 构建 的 CE - 318 大 气 水 汽 改 
进 型 反 演 模型 [公式 (5) ] 。 经 模型 参数 的 相关 性 分 
析 及 精准 度 验 证 ,3 个 模型 的 反 演 结果 都 较为 理想 。 
其 中 ,基于 CE -318 的 大 气 水 汽 反 演 模型 的 反 演 结 
果 各 项 检验 指标 均 为 最 佳 , 故 选 取 其 为 研究 区 域 的 


大 气 水 汽 是 反映 大 气 环境 变化 的 一 个 重要 因 


最 优 模型 ,表明 该 模型 可 为 石家庄 地 区 大 气 水 汽 的 
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反 演 监测 提供 新 的 思路 和 方法 。 

本 文 建 模 数据 为 点 状 数据 源 时 刻 , 故 基 于 单 站 
CE -318 实测 数据 ,并 加 入 探 空 这 一 变量 修正 参数 
构建 的 改进 型 模型 所 得 结果 为 点 状 大 气 水 汽 值 。 可 
依照 降水 前 大 气 水 汽 含量 会 有 快速 升 高 过 程 的 现 
象 ,根据 点 状 水 汽 值 制作 该 站 点 的 大 气 水 汽 时 间 
序列 图 ,用 来 进行 水 汽 日 变化 .季节 变化 及 年 际 变 
化 等 时 间 方 向 上 的 气候 分 析 。 同 时 ,也 可 根据 需 
求 利 用 合适 的 空间 插值 法 将 点 状 大 气 水 汽 数据 拟 
合 为 面 状 结果 ,制作 区 域 水 汽 分 布 图 ,有 利于 水 汽 
分 布 的 动态 变化 检测 ,以 揭示 天 气 状况 形成 发 展 
和 移动 的 规律 。 

由 于 本 文 在 研究 时 选取 UTC 0 时 作为 3 种 数据 
的 可 对 比 时 刻 , 这 必然 会 造成 其 他 时 间 有 效 数据 的 
大 量 牺牲 ,对 模型 的 建立 和 使 用 增加 了 时 间 上 的 限 
制 条 件 。 在 之 后 的 研究 中 ,可 通过 整理 历史 数据 以 
增加 其 他 时 刻 或 改进 数据 筛选 与 处 理 的 方法 ,最 大 
限度 地 提高 所 获 数据 的 利用 率 ,以 期 提高 模型 的 精 
度 和 拓展 模型 的 普 适 性 。 
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Inversion Model for Atmospheric Water Vapor over the Shijiazhuang Area 
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(1. Hebei Province Laboratory of Environmental Evolution and Ecological Construction , Shijiazhuang 050024 , Hebei , China ; 
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Abstract: It is of important practical significance for effectively assessing the functions of atmospheric ecological 
services and strengthening the construction of ecological environment to study the inversion of atmospheric water va- 
por and understand its spatiotemporal change. Based on the CE —318 observed atmospheric vapor data , sounding 
data and surface vapor pressure data ,the traditional regression model, CE — 318-based atmospheric water vapor in- 
version model and seasonal model were constructed by traditional regression analysis and improved regression analy- 
sis. The optimal inversion model of atmospheric vapor over the study area was selected by the accuracy evaluation 
analysis. The main conclusions are as follows: The inversion results of the three models constructed were all rela- 
tively satisfactory , and the inversion results of the atmospheric water vapor inversion model based on CE —318 were 
the most accurate. Its item inspection indexes were the best values in the accuracy tests of the three models , and the 
atmospheric water vapor inversion models were optimal for the study area. 
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